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First Passage Percolation and Competition Models 

Nathaniel D. Blair-Stahn 
April 20, 2007 

Abstract 

This paper is a survey of various results and techniques in first passage percolation, 
a random process modeling a spreading fluid on an infinite graph. The latter half of 
the paper focuses on the connection between first passage percolation and a certain 
class of stochastic growth and competition models. 
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1 Introduction 

First passage percolation is a random process on a (typically infinite) graph. Hammersley 
and Welsh JHW65] introduced first passage percolation as a model of fluid flow through a 
randomly porous medium. In this model, each edge e in the graph is assigned a random 
nonnegative number T(e), called the passage time of e, which is interpreted as the time it 
takes to cross the edge in either direction. In other contexts, r(e) may represent a weight or 
a capacity, but we shall stick with the passage time interpretation. The picture to keep in 
mind is that of a fluid emanating from some source vertex and flowing outward through the 
edges of the graph according to the prescribed passage times. Equivalently, one can think 
of an infection spreading out from some initial locus and transmitted between neighboring 
vertices at random times. 

First passage percolation can be defined on any graph, but the most commonly studied 
model is the one in which the underlying graph is the integer lattice TL^ for d>2. This is 
the model we will focus on, although we briefly discuss models on certain random infinite 
graphs in § 12.61 The most basic results in first passage percolation rely on ergodic theory 
and the theory of subadditive processes. In fact, the study of first passage percolation was 
an impetus for the development of Kingman's subadditive ergodic theorem [ Kin68j . A good 
general reference detailing the fundamental results in first passage percolation is |Kes86] . 

Based on heuristic arguments, the growing interface described by first passage perco- 
lation is believed to belong to the Kardar-Parisi-Zhang (KPZ) universality class |KPZ86j 
of models in statistical physics. In particular, the Eden growth model |Ede61] . which can 
be described in terms of a particular first passage percolation process, has been studied 
extensively in this context. We will define the Eden growth model in § 13.11 In § 12. 41 we will 
see some of the progress that has been made in verifying various predictions from statistical 
physics. 

The rest of the paper is organized as follows. In § [2], after giving the precise definition 
of the first passage percolation process and introducing some of the topics of interest, we 
summarize the early results in the subject and proceed to describe some of the more recent 
work that has been carried out. In § [3] we describe how to define growth processes and 
competition models based on first passage percolation and summarize recent work in this 
direction. 



2 First passage percolation 

2.1 Basic definitions 

Let U^ be the integer lattice of dimension d > 2, which we consider both as a graph and 
as a subset of M'^. Two vertices u,v E Z'^ are adjacent if ||m — v\\i = 1, and we denote the 
edge set of Z'^ by E{l/). Let {T(e)}eg£;(zd) be a collection of nonnegative random variables 
indexed by the edges. We call r(e) the passage time of the edge e, and it represents the 
time needed to cross the edge in either direction. The joint law of the passage times r(e) 
determines the first passage percolation process. 

In order to do anything useful with the first passage percolation model, we need to 
make some assumptions about the distribution of the passage times. Typically, the minimal 
assumption one makes is that the passage times {T{e)}e^E('Ld) are stationary and ergodic 
with respect to translations of 7/. More explicitly, we can consider the canonical sample 
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space ri = (]R+)^'^'^"\ equipped with some probability measure v defined on the product 
cr- field. If w G f2 is a realization of z/, the passage times for w are given by T(^(e) = w(e). 
Now, for each m G Z'^, let 6'^ : fi — )■ fi be the natural shift operator defined by 

6'„a;(e) = a;(e + u), 

where the notation "e + w" has the obvious meaning. The passage times r(e) are stationary 
if for each u, the measure u is ^^-invariant (i.e. z/ o 9~^ = z/). Stationary passage times 
are ergodic if any event which is invariant under every 6u has probability or 1 (i.e. any 
event A G Q such that 0^^^ = ^ for all u must satisfy z/(y4) G {0, 1}). Throughout the 
rest of the paper we will assume that the passage times are stationary and ergodic, and 
we will be most interested in the case where they are in fact independent and identically 
distributed (i.i.d.) and have finite expectation. In particular, our focus in § [3] will be on 
i.i.d. exponential passage times. 

We now define several concepts that will be discussed in more detail in later sections. 
Suppose that {T(e)}eg£;(zd) is a collection of passage times describing a first-passage perco- 
lation process on Z"^. If 7 is a path in Z'', then the passage time of 7 is 

T(7) = 5^r(e). 

li U,V C Z*^, the passage time from f/ to V is 

T(f/, V) = inf {T(7) : 7 is a path from U to V}. (2.1.1) 

(If t/ or y is a singleton, we will write its unique element in place of the set when using 
this notation or other similar notation.) We can extend this definition to subsets of M'' as 
follows: If A C M*^, let A consist of all the lattice points that are closest to A, i.e. 

A = Iv e Z'^ : V e X + [-^,lY ^°^ ^^"^^ X G a| , 

and for U,V C M^ set T(f/, V) := T{U, V). 
For example, for each n = 0, 1, 2, . . ., let 

ri = (ra,0, ...,0) G Z'^ and Hr, = {z e Z"^ : Zi = n}. 

We refer to T{0,n) as a point-to-point passage time and T(0, i^„) as a point-to- 
hyperplane (or point-to-line when d = 2) passage time. We will see in § 12.21 that the 
passage times T(0, n) and T(0, iJ„) satisfy a law of large numbers (Theorem 12.21) . which 
shows that first passage percolation has an asymptotic speed along the coordinate axes. 

One of the primary objects of interest in first passage percolation is the set B{t) of 
vertices that can be reached from the origin by time t, or a continuum version B{t) of this 
set in which each v G B(t) is replaced with a unit cube centered at v. That is, 

B{t) = {v eZ'^ : T(0, v) < t} and B{t) = {x eM.'^ : T{6, x) < t}. 

If we think of the percolation process as modeling an infection spreading outward from 
the origin, B{t) is the set of vertices which are infected at time t. (We will return more 
explicitly to this interpretation of B{t) in § [31 where we discuss growth and competition 



models.) One of the fundamental results about first passage percolation is that, under some 
mild hypotheses for the passage times, B{t)/t converges almost surely to a deterministic 
shape (Theorem I2.4p . so the process in fact has an asymptotic speed in all directions 
simultaneously. We will discuss this so called shape theorem further in § 12. 3[ and in § 12.41 
we will discuss the related question of deviation bounds for the convergence. 

Another topic that arises naturally in the study of first passage percolation is that of 
time-minimizing paths, or geodesies. For sets U,V G M"', if T{U, V) = T{'~f) for some (nec- 
essarily finite) lattice path 7 (i.e. 7 achieves the infimum in (I2.1.ip ). we call 7 a geodesic 
from U to V, and we denote any such path by G{U, V). More generally, a finite or infinite 
path 7 in Z'^ is called a geodesic if every finite subpath 7' of 7 satisfies ^(7') = T{u', v'), 
where u' and v' are the endpoints of 7'. Observe that a finite path 7 with endpoints u and 
u is a geodesic if and only if 7 = G{u, v). When the passage times are i.i.d., it is easy to see 
that G{u,v) exists and is unique a.s. for each u,v gU^ if and only if r(e) is a continuous 
random variable (see e.g. |WW98t Lemma 8]). See |GM05t § 4] for conditions guaranteeing 
the existence and uniqueness of finite geodesies with stationary passage times. We will 
discuss the existence of infinite geodesies in § 12. 5[ 

2.2 The subadditive ergodic theorem and the time constant 

One property of the point-to-point passage times {T(m,n)}o<m<n that is immediate from 
the definition (12.1.11) is 

T(0, n) < T(0, m) + Tim, n) for all < m < n. 

This motivates the following definition: A doubly indexed process {Xm,n}o<m<n is called 
subadditive if Xo,„ < Xo^m + X^^n for all < m < n. 

The main result about subadditive processes is the subadditive ergodic theorem, which 
was developed by Kingman |Kin68j to study point-to-point passage times and is now a 
standard tool in first passage percolation and other applications. The following version. 



due to Liggett Lig85a| , is an improvement on Kingman's original result. (Instead of (b) and 
(c) below, Kingman assumed that the distribution of {Xm+k,n+k\ Q<:_rn<n ^'^^^ i^o^ depend 
on fc, in which case it follows from (a) that X^^n < X^ „ -|- Xm,n for aX\ C, <m < n.) 



Theorem 2.1 (Subadditive ergodic theorem |Kin68j . Lig85a] ). Suppose X^.n, < m < n, 

is a family of random variables satisfying 

(a) Xo,„ < Xo,m + Xm,n for all < 171 < 71. 

(b) For each k > 1, the sequence {Xnk,{n+i)k}^^Q is stationary. 

(c) The distribution of the sequence {Xm,m+k}k>i does not depend on m. 

(d) EXo+i <oo. 
Then 

(i) lim„_j.oo EXo,„/n = inf„ EX^n/n = 7 for some 7 > —00. 

(ii) The limit X = lim„^oo Xo,ri/n exists and is less than +00 a.s. 

(iii) // there is some c < 00 such that E X^^ > —en for all n, then the convergence in (ii) 
also holds in L} , so EX = 7. 



(iv) If the stationary sequences in (b) are ergodic, then X = 7 a.s. 

If the random variables in Theorem 12. II are all degenerate, then the theorem reduces to a 
well-known result about subadditive functions (see e.g. |Kal02t p. 191] or |Kin73] ) . On the 
other hand, if {^fc}fc>i is a stationary sequence with E |,^fc| < 00, then Xm,n := ^m+i + - • ■+C,n 
satisfies (a)-(d) and (iii) (with equality in (a), so the process {Xm,n} is in fact additive), 
and Theorem 12.11 reduces to Birkhoff's ergodic theorem in this case. 

If the passage times r(e) are stationary and ergodic with finite expectation, the random 
variables Xm,n '■= T{rh,n), < m < n, are easily seen to satisfy (a)-(d), (iii), and (iv) of 
Theorem 12. H so T{0,n)/n converges a.s. to some constant yUi < 00. The constant ^i is 
known as the time constant in the direction 1, and its reciprocal is the the asymptotic 
speed of the first passage percolation process along the coordinate axes. When the passage 
times are i.i.d., it turns out that the scaled point-to-hyperplane passage times T(0, if„)/n 
converge to the same constant fii. This was first proved by Wierman and Reh |WR78] . 
and can be deduced from the shape theorem (Theorem 12.41 below - see |Kes86t pp. 166- 
167]). (Note that the process T^rh, Hn) is not subadditive, so we cannot apply Theorem 12. II 
directly.) We restate these two results for i.i.d. passage times in the following theorem. 

Theorem 2.2 (Time constant |Kin68j . |WR78j ) . Suppose the passage times {T{e)}e^x'i are 
i.i.d. with finite expectation. Then there is a constant ;Ui < cxd such that 

Finally, we mention two basic results about the time constant /ii in the case of i.i.d. 
passage times. First, it is easy to see that in general /ii < E[r(e)] (see |HW65l Theo- 
rem 4.1.9]). Also, observe that /ii = corresponds to infinite percolation speed, so that the 
process has superlinear growth. There is a simple criterion for deciding when this occurs. 

Proposition 2.3. For i.i.d. passage times T{e), the time constant ^i is nonzero if and only 
z/Pr[r(e) = 0] < pc{1j'^), where Pc(^'^) is the critical value for Bernoulli bond percolation 
onU^. 

A proof of Proposition 12.31 can be found in |Kes86l § 6]. A heuristic argument goes 
as follows. If Pr[r(e) = 0] > pJ(L^\ then there is a.s. an infinite cluster in Z"^ on which 
the travel time between any two vertices is zero. It will a.s. take only finite time to reach 
this cluster from the origin, at which point the process can head off in any direction with 
infinite speed. On the other hand, if Pr[r(e) = 0] < pJ(L'^\ then a.s. all the clusters on 
which infinite speed can occur have finite size. Thus, the process can only travel a finite 
distance before it has to step off one of these clusters and accumulate some positive travel 
time before reaching the next cluster. It is not too hard to show that this accumulated 
travel time must with high probability increase linearly with the distance traveled, so that 
the asymptotic speed is finite a.s. The situation at the critical value pJ(L'^) is a bit more 
delicate, but Proposition 12.31 shows that the asymptotic speed in this case is infinite. 

2.3 The shape theorem 

It is natural to generalize the idea of the time constant and consider the speed of percolation 
in arbitrary directions rather than just along the coordinate axis. In particular, for any 
X G M'^ with rational coordinates, we can apply Theorem 12.11 to see that there is some 



constant /i(x) such that T{0,nx)/n — )■ fi{x) a.s. With this notation we have /ii = fi{l). 
For i.i.d. passage times, it is not difficult to show that the function /x : Q"^ — )■ [0, oo) 
is Lipschitz continuous and hence can be extended to all of W^, and that the resulting 
function /i : M"' — )■ [0, oo) is either identically zero or defines a norm on M°'. We will refer to 
fi as the norm for the ffist passage percolation process when appropriate; more generally, 
we will refer to fi as the shape function for the process because of its role in the shape 
theorem, which we now describe. 

Recall the definitions of the growing shapes 

B{t) = {v eZ"^ : T(0, v) < t} and B{t) = {x eW^ : T(0, x) < t}. 

Under some moment conditions on the passage times, it can be shown that B{t)/t converges 
almost surely to the unit /i-ball So = {a; G M'^ : /i(x) < 1} as t — )■ oo. This result is known 
as the shape theorem, and an in probability version was first proved by Richardson [Ric73] 
for d = 2. Cox and Durrett |CD81j used a result of Kesten (found in |Kin73l p. 903]) to 
strengthen Richardson's result to an almost sure version. The following version, valid in 
any dimension, is proved by Kesten in |Kes86] . 

Theorem 2.4 (Shape theorem |Kes86t Thm. 1.7]). Suppose that {T{e)}f,^E{z^) (^''^^ i.i.d. 
passage times such that Emin{r(ei)'^, . . . , T{e2dY} < oo (where {ei, . . . , 62^} is any set of 
2d distinct edges). Let fi be the shape function for the process, let /ii = fi{l) be the time 
constant, and let B^ = {x E M.'^ : fi{x) < 1}. 

1. If fii > 0, then Bq is compact and convex with nonempty interior, and for any e > 0, 

(l-6)i?oC^C(l + 6)i?o 

for all large t almost surely. 

2. If fii = 0, then /i = (so Bq = M^^, and for any compact set K C W^, 

t 
for all large t almost surely. 

The moment condition on the passage times in Theorem 12.41 is optimal, in the sense 
that if it fails then limsup^^ooT(0, v)/||v||i = 00 a.s. However, even without any moment 
conditions on the passage times r(e), it is possible to define modified passage times T(m, v) 
for u,v eZ/ and a corresponding set B{t) such that an analogue of Theorem 12.41 holds (see 
|Kes86]l By Proposition [Ml we see that Bq = M'^ if and only if Pr[r(e) = 0] < pdZ"^). The 
convexity of Bq follows from subadditivity, and when Bq 7^ M"', compactness and nonempty 
interior follow from the fact that /i is a norm. Otherwise, little is known about the limit 
shape Bq other than the obvious fact that it must have all the symmetries of Z"'. Kesten 
|Kes86l § 8] shows that if the passage times are i.i.d. exponential and d is large, then Bq 
is not a Euclidean ball, casting doubt on the conjecture that Bq might be a disc for d = 2 
based on early Monte Carlo simulations [EdeGlj . Durrett and Liggett |DL81j show that 
there are i.i.d. passage times for which Bq has flat edges but is not a diamond or a square. 
In particular, this occurs if r(e) is nontrivial but attains some nonzero minimum value with 



probability greater than p^"'(Z'^), where p^"'(Z'^) is the critical value for directed Bernoulli 
bond percolation on Z'^. 

There is also a version of the shape theorem for stationary passage times. Boivin |Boi90] 
proves that if the passage times r(e) are stationary, ergodic, and have finite moment of 
order c? + e for some e > 0, then B{t)/t converges a.s. to a deterministic shape Bq. In the 
stationary case, the shape function fi may take on both zero and strictly positive values so 
that the limit shape Bq can be an unbounded proper subset of M'^. However, if /i(x) > for 
every unit vector x, then Bq is compact, convex, has nonempty interior, and is symmetric 
with respect to reflection through the origin. {Bo may fail to have further symmetries since 
isotropy may not hold in the non-i.i.d. case.) Conversely, Haggstrom and Meester |HM95] 
show that any set Bq C M'^ with these properties can arise as the limit shape for some 
collection of stationary passage times. 

2.4 Deviations in the passage times and the growing shape 

Throughout this section we will assume that the passage times {T{e)}^(zE(z'') are i.i.d. and 
satisfy the hypotheses of Theorem 12.41 so that B{t)/t — )■ Bq a.s. We further assume that 
Pr[r(e) = 0] < Pc{'^'^) so that /ii > and the limit shape Bq is compact. 

2.4.1 The variance of T(0, n) 

Once we know that B(t) converges, we can ask how much it deviates from the limit shape 
Bq. There are various ways to approach this problem. As a first step, we consider the 
variance of T(0, n). It is predicted that the standard deviation of T(0, n) is of order n^ 
for some constant x = x{d)- Based on heuristic arguments from statistical physics, it 
is expected that x(2) = 1/3 (see e.g. |KS91j . [KPZ86J ). This conjecture is supported 
by simulations and by rigorous results for related growth models (e.g. |BDJ99j . |JohOOa] . 
[JohOObj ). which we shall discuss in § 12.71 The situation is less clear for higher dimensions 
d, although it is generally believed that x is nonincreasing in d (see |NP95] for a discussion). 
So far, the only general bound on x, due to Kesten, is xi.^) < 1/2 for all d: 

Theorem 2.5 (Kesten |Kes93] ). //E[r(e)^] < oo, then there are positive constants ci and 
02 such that 

Ci < Var[T(0, n)] < C2n. 

Kesten proves Theorem 12.51 using martingale methods (the "method of bounded differ- 
ences"). Although Theorem 12.51 provides the best known bounds for a general distribution 
on the passage times, better bounds have been proved for certain classes of distributions. 
For example, Benjamini, Kalai, and Schramm |BKS03j use an inequality of Talagrand 
|Tal94| Thm. 1.5] to show that Var[T(0, n)] = 0(n/logn) if the passage times have the 
uniform distribution on {a,b}, where < a < 6. They note that the essential feature of 
first passage percolation needed to prove both their result and Kesten's is that the number 
of edges e G E{l/-) such that modifying r(e) increases T(0,n) is bounded by a constant 
times n. Building on the methods in |BKS03j , Benaim and Rossignol |BR06aj use a Gaus- 
sian version of Talagrand's |Tal94] inequality and apply the techniques of |BKS03] to prove 
0{n/ logn) variance for a large class of i.i.d. absolutely continuous passage times, including 
exponential. 

As for lower bounds on the variance, Newman and Piza [NP95J prove that in dimen- 
sion d = 2, Var[T(0,'n)] = Q{logn) under certain hypotheses on the passage times. In 
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particular, if we set A = mi{x : Pr(r(e) < x) > 0}, the condition assumed in [NP95J is 
that 

Pr(r(e) = A)<p(A), (2.4.1) 

where 

>,(Z2) if A = 0, 

pf (Z2) if A > 0. 



PW = ^ ^dir,"y2^ 



Based on Proposition 12.31 and the results in [DL81J . this condition is necessary for the 
shape Bq to be compact and for its boundary to have no fiat edges; it is suspected that 
(12.4.11) should also be sufficient for this to hold (see |NP95] ). Pemantle and Peres |PP94] 
use different methods to prove ^^(logn) variance for the special case of exponential passage 
times in d = 2. On the other hand, for higher dimensions d, it is still not known whether 
the variance of T(0, n) even diverges as n — > oo. 

2.4.2 Large deviation bounds for T(0, n) 

If the assumption of finite variance for the passage times is strengthened to the existence of 
a finite exponential moment, then one can obtain good bounds on the deviation of T(0, n) 
from its expected value, and on the deviations of E[T(0, n)] from nfi. The following theorem 
is due primarily to Kesten |Kes93] , with the upper bound in (12.4.31) being an improvement 
made by Alexander. 

Theorem 2.6 (Kesten [Kes93j . Alexander |Ale93 j). //E[e'^'^'^^^] < oo for some 7 > 0, then 
there exist positive constants 01,02,03,04,05, such that 



Pr 

and 



T(0,n)-E[T(0,n)] 



n 



> X < 010-^2"^ for X < Can, (2.4.2) 



04- < E[T(0, n)] -nfi< c^n^^^ log n. (2.4.3) 

n 



Note that the lower bound in f l2.4.3p strengthens the trivial inequality E[T(0,n)] > nfi 
implied by Theorem 12.11 Both Theorems 12.51 and Theorem 12.61 remain valid if T(0, n) 
is replaced by T(0,H„) (see |Kes93j or |Ale93j ). or if n is replaced by any t> G Z"' and 
n is replaced by \\v\\i. In fact, using versions of (12.4.21) and (12.4.31) valid for arbitrary 
directions, Kesten |Kes93t Theorem 2] shows that there is some constant C (depending on 
the dimension d and the distribution of T(e)) such that almost surely, 

1-^-^j j.BoC^C^l + -^j.So foralllarget. 

Some improvements of Theorem 12.61 are available in certain situations. Talagrand |Tal95t 
§8.3] shows that the upper bound in (I2.4.2p can be strengthened to 0{e~'^^ ) if E[T(0, n)] 
is replaced with a median of T{0,n). For the same class of distributions considered in 
[BR06a] (with the added assumption of finite exponential moment), Benaim and Rossignol 
[BR06bJ prove that (I2.4.2p still holds if the y/n in the denominator is replaced by -y/n/logn. 
Instead of the Talagrand- type inequalities used in |BR06a] , the techniques used in |BR06b] 
involve modified Poincare inequalities arising from the context of "threshold phenomena" 
for Boolean functions. 



2.4.3 Scaling exponents for the growth process 

In the statistical physics literature (see e.g. |KS91] ). the fluctuations of a randomly growing 
shape such as B{t) are studied in terms of two exponents x ^i-nd ^, which describe respec- 
tively the longitudinal and transverse fluctuations in the surface of B{t). For example, it 
is expected that the standard deviation of the time T(0, H) at which B{t) flrst reaches a 
hyperplane H at distance r from the origin is of order r^, while the set of points in H 
which are likely to be flrst reached by B(t) is expected to have diameter on the order of r^. 
There are various ways to deflne x ^i^nd ^ precisely, and it is an open problem to determine 
whether the various deflnitions are equivalent. 

The exponents x ^i-nd ^ are not expected to depend on the underlying distribution of 
the T(e)'s, at least under certain hypotheses (for example, 02.4. ip above - see |NP95j or 
[LNP96J ). A priori, x cind ^ could depend on the direction of travel, but it is expected 
that they should be the same in any direction in which the boundary of Bq has nonzero 
curvature, at least in low dimensions. The values of x and ^ are expected to depend on the 
dimension d, but heuristic arguments suggest that the scaling identity x = 2.^ — 1 holds in 
all dimensions (see |KS91j ). As noted in the introduction, flrst passage percolation models 
are expected to belong to the KPZ universality class |KPZ86j . leading to the prediction 
that x(2) = 1/3 and (in accordance with the scaling identity) ^(2) = 2/3. 

We now describe some of the progress that has been made towards computing the 
exponents x ^^^ C Since x ^^^ ^ might depend on the direction of travel, we will write 
Xx and C,x to denote their values in the direction of some unit vector x G M"^. In |NP95j . 
Newman and Piza show that in any dimension d, if (12.4.11) holds, then Xx > (l~('^~l)^x)/2 
(this was proved by Wehr and Aizenman |WA90j for d = 2). Then they show that, under the 
same hypothesis (I2.4.ip . if the passage times have flnite exponential moment, then C,x < 3/4 
for any x which is a direction of curvature for Bo (i.e. a direction in which the boundary of 
Bq has nonzero curvature). For d = 2, this yields Xx > 1/8 in any direction of curvature x, 
improving the previously mentioned logarithmic lower bound on Var[T(0, nx)]. It is easy 
to show that any compact convex set has a direction of curvature |NP95t Lemma 5], so in 
d = 2 there is at least one direction x such that Var[T(0, nx)] = n(n^^^) when the r(e)'s 
have flnite exponential moment. 

The method in |NP95j used to prove C,x < 3/4 makes use of an exponent x' analogous 
to Xy but which also takes into account the deviations of E[T(0, nx)] from nfi{x). The 
Kesten-Alexander deviation bounds (Theorem 12.60 imply that x' ^ 1/2- Newman and 
Piza then use a rigorized version of the heuristic argument from |KS91j for the scaling 
identity x = 2i^ — 1 to show that ^£ < (1 + xlf^y which yields the bound ^x < 3/4. 

In |LNP96j . Licea, Newman, and Piza extend the methods in |NP95j to obtain lower 
bounds on various versions of the exponent ^. Combining the trivial bound x ^ with the 
(nonrigorous) scaling identity yields the nontrivial bound C, > 1/2, which is expected to 
hold in all dimensions. The value ^ = 1/2 corresponds to what is called a diffusive process, 
and it is believed that, at least in low dimensions, flrst passage percolation should in fact 
be superdiffusive, i.e. ^ > 1/2. Using progressively weaker deflnitions ^'•^\ ^'•^\ ^'•^^ for ^, 
Licea, Newman, and Piza prove 

e\d)>l/{d+l), e^\d)>l/2, and e^\2)>?,/h, 

assuming that the passage times satisfy (12.4.10 and/or E[r(e)^] < oo. The latter two 
bounds correspond to superdiffusivity as predicted by the physical models. While the flrst 
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bound is subdiffusive, it is nontrivial from a mathematical perspective, and may be useful 
because the exponent ^'•^^ has certain advantages over the other two definitions of ^. 

2.5 Infinite geodesies 

Recall that a geodesic is a time-minimizing path in first passage percolation, and that 
G{U, V) denotes a geodesic between the sets U and V when such a path exists. Suppose 
that G(m, v) exists and is unique for each pair of vertices u,v & Z'^. For any u G Z'^, we 
define the tree of infection of u, ^{u), to be the (graph theoretic) union of all the finite 
geodesies starting at u: 

T{u)= \jGiu,v). 

The fact that T{u) is a tree follows from the uniqueness of the geodesies. If we think of the 
percolation process as modeling an infection spreading outward from m, then the unique 
path in T{u) from u to another vertex v traces the route by which v became infected. 

Let K{T{u)) denote the number of topological ends in T{u) - that is, the number of 
semi-infinite paths in T{u) starting at u. We call any such path a one-sided geodesic 
starting at u. A standard compactness argument shows that K{T{u)) > 1 for any u. (The 
set of finite geodesies starting at u can be viewed in a natural way as a compact space, so 
it must contain a limit point since it has infinitely many elements.) In [New95j . Newman 
uses the Kesten- Alexander deviation bounds (Theorem 12.61) and methods similar to those 
in |NP95j to show that if the passage times are i.i.d. and the curvature of the boundary 
of Bo is uniformly bounded away from 0, then K{r{u)) = oo a.s. for any u. While the 
assumption of uniform curvature is plausible, there are no i.i.d. probability measures on the 
passage times for which Bq is known to have this property. Hoffman |Hof05bj has shown 
that K{r{u)) = oo under a much weaker assumption on the limit shape, namely that Bq is 
not a polygon. Although there are no i.i.d. passage times which are known to satisfy this 
assumption either, the result of |HM95j shows that there are stationary passage times for 
which it holds. We will revisit this topic in § 13.31 where we discuss the connection between 
the existence of geodesies and the question of coexistence in a certain competition model 
obtained as a projection of first passage percolation. 

Newman proves the above result by considering geodesies with an asymptotic direction. 
If X is a unit vector in M"' and 7 is a one-sided geodesic with vertices fo,fi,f2, . • ., then 
7 has asymptotic direction x if lim„_j.oo'yn/||'yn||2 = x, and we call 7 an a;-geodesic. It is 
not known in general whether x-geodesics exist or whether every one-sided geodesic must 
have a direction, but Newman |New95j gives affirmative answers to both questions under 
the assumption that Bq is uniformly curved. Licea and Newman |LN96] use some of the 
ideas in |New95j to prove a uniqueness result for x-geodesics when d = 2: 

Theorem 2.7 (Licea and Newman |LN96] ). Suppose the passage times {r(e)}egz2 are i.i.d. 
with continuous distribution. Then for Lebesgue almost every x on the unit circle, 

Pr (There exist disjoint x-geodesics) = 0, 

and hence any two x-geodesics must coalesce. 

Naturally, in addition to one-sided geodesies, one can also consider two-sided geodesies, 
i.e. geodesies which are infinite in both directions instead of just one. In contrast to one- 
sided geodesies, it is unclear whether two-sided geodesies exist at all. (Since a two-sided 
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geodesic does not have a fixed starting point, the compactness argument used to prove the 
existence of one-sided geodesies fails in this case.) Most of the work on two sided-geodesics 
has focused on the case of continuous i.i.d. passage times in d = 2, where it is expected 
that two-sided geodesies do not exist. We summarize some results in this direction. 

In jWeh97] . Wehr shows that almost surely, the number of 2-sided geodesies is either 
or oo in (i = 2, and that an analogous result holds in d dimensions for locally weight- 
minimizing hypersurfaces instead of curves. This result is equivalent to the statement that 
the number of ground states in the random exchange Ising model (REIM) is 2 or oo a.s. 

Using Theorem 12. 7[ Licea and Newman [LN96J show that for Lebesge-a.e. unit vector 
X G M?, there cannot exist an (x, — i:)-geodesic, i.e. a two-sided geodesic with asymptotic 
directions x and —x. Wehr and Woo |WW98j show that if H is any half-plane in M^, then 
there can exist no two-sided geodesies contained entirely within H. As a corollary, any 
two-sided geodesic must intersect every line with rational slope. 

2.6 Isotropic models of first passage percolation 

As we have seen on various occasions above, one disadvantage of the first passage perco- 
lation model on Z'^ is that we do not have much information about the limit shape Bq. 
One way to get around this is to define a stochastically isotropic model so that symmetry 
considerations imply that Bq must be a Euclidean ball. For example, Vahidi-Asl and Wier- 
man |VAW90j introduce models in which the underlying graph is either a random Voronoi 
tesselation of the plane or its dual Delaunay triangulation, where the centers of the Voronoi 
cells are given by a Poisson point process on M^. Howard and Newman |HN97j introduce 
a different model, in which the underlying graph is the complete graph with vertices given 
by a Poisson point process on M'' and the passage times are given by r(e) = |e|", where 
a > 1 and |e| denotes the Euclidean distance between the endpoints of e. 

Using a random graph for the for the percolation process introduces various technical 
problems, but nevertheless, versions of many of the results familiar from the Z'^ model 
still hold for these Euclidean models. For example, in the Voronoi and Delaunay models, 
there is a time constant |VAW90] . a shape theorem |VAW92] . and deviation bounds similar 
to those in Theorems 12.51 and 12.61 |Pim05] . Furthermore, since Bq has uniform curvature 
in this model, Pimentel |Pim04j is able to use the techniques in |NP95j and |New95j to 
show that the transversal fiuctuation exponent ^ < 3/4, and that almost surely, every one- 
sided geodesic has an asymptotic direction and there exists a one-sided geodesic in every 
direction. Similar results are proved by Howard and Newman for their model in |HN97] . 
|HN99] . and |HN01] . We also mention that jPim04] contains results about a competion 
model on the Delaunay triangulation analogous to the competition model on Z'^ described 
in § 13.31 below. 

2.7 Directed first passage and last passage percolation 

The process we have been referring to as first passage percolation is more properly called 
undirected first passage percolation. One can also consider directed first passage percolation 
or a related model called (directed) last passage percolation. Both of these models are 
defined similarly to undirected first passage percolation, except that only increasing paths 
(defined below) are allowed. Certain versions of directed last passage percolation are much 
better understood than undirected first passage percolation. We now describe the directed 
models and summarize some of the most interesting results. 
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Let {T(e)}egZ'^ be a collection of i.i.d. nonnegative passage times. A path 7 in Z'^ is 
called increasing if each step in 7 is made by increasing a single coordinate by 1. For 
M, w G Z'^, write u < v ii Ui < Vi for 1 < i < d. If m < f , define the (directed) first-passage 
time from m to f to be 

Tminiu,v) = min{T(7) : 7 is an increasing path from u to v}, 

and define the last-passage time from n to f to be 

Traax{u,v) = max{T(7) : 7 is an increasing path from u to v}, 

where T(7) is defined as in the undirected case. As before, we can extend these definitions 
to passage times between two points in M*^. The directed models are often defined with 
passage times t{v) on the vertices v of Z'^ rather than the edges, but the analysis is similar 
with either convention, so we will stick with edge passage times. 

As with the undirected first-passage times, the directed first-passage times are subad- 
ditive, whereas the last-passage times are superadditive, i.e. for vertices u < v < w we 
have 

Applying Theorem 12. II to the first-passage times and a superadditive version of Theorem 12. II 

to the last-passage times implies that there are shape functions g, h : (M_|_)'^ — )■ [0, cxd) such 
that, for all x G (M+)'^, 

lim = g[x) a.s. and lim = h[x) a.s. 

n— >co n n— >-oo n 

Furthermore, we can define growing shapes analogous to B{t): 

U{t) = {xe (M+)^ : r„,in(0, x) < t} and V{t) = {x E {^+Y : T^ax(0, x) < t}. 

Under appropriate conditions on the distribution of T(e), Martin |Mar03j proves a shape 
theorem for the directed models: 

U{t)/t -^ Uo a.s. and V{t)/t -^ Vq a.s., 

where Uq = {x : g{x) < 1} and Vq = {x : h{x) < 1}. In the first-passage case, subadditivity 
implies that Uq is convex, whereas in the last-passage case, superadditivity implies that 
(M+)'^ \ Vq is convex. 

In contrast with the undirected model, there are two special cases of directed last- 
passage percolation in d = 2 for which the shape function h{x) is known explicitly. A 
theorem of Rost [RosSlj (see also |BS04] . |Mar03j ) implies that for exponential passage 
times with mean 1, 

h{x) = ||x||i/2 = iy/x^+ a/x^)^- 

Johansson |JohOOa] shows that for geometric passage times with parameter q, 

, . X , . X q{xi+X2) + 2^X1X2 

h[x) = hq[x) = . 

1-q 

These are the only two nontrivial cases where the shape function for i.i.d. passage times 
is known, in any of the directed or undirected, first- or last-passage models. However, 
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Seppalainen |Sep98] finds the limiting shape for a particularly simple stationary model of 
directed first-passage percolation on (Z_(_)^, in which vertical edges have a deterministic, 
constant passage time and horizontal edges have i.i.d. Bernoulli passage times. 

In fact, Johansson |JohOOa] not only identifies the shape function in the i.i.d. geometric 
last-passage model, but extends the techniques of Baik, Deift, and Johansson [BDJ99J to 
show that the passage times T^sx{0,nx), appropriately centered and scaled, converge in 
distribution to the Tracy- Widom [TW94J distribution for the largest eigenvalue in a random 
matrix sampled from the Gaussian Unitary Ensemble (GUE). In particular, it is shown that 
the standard deviation of Tmax(0, nx) in this model is of order n^/^, so that x = 1/3 in 
accordance with the predictions of KPZ universality |KPZ86| . 

Finally, we mention another model which can be viewed as a continuum version of the 
directed last-passage percolation model defined above, and in fact can be obtained as a 
limit of last-passage percolation with i.i.d. geometric passage times (see, e.g. |Joh02] ) . This 
model was introduced by Hammersley |Ham72] as a method for approaching Ulam's prob- 
lem |Ula61] of finding the distribution of the longest increasing subsequence in a random 
permutation. 

Consider a unit-rate Poisson process on M^. Analogous to increasing lattice paths, we 
can define an increasing path between Poisson points to be a path 7 that moves only up and 
to the right. That is, 7 is a sequence of Poisson points such that if x and x' are consecutive 
points in 7, then x < x' . We define the length of 7 to be the number of Poisson points 
it contains. Let L(r) be the length of the longest increasing path between Poisson points 
contained in the square [0,r]^. Conditional on the event that [0,r]^ contains N points, 
Lir) has the same distribution as the length of the longest increasing subsequence in a 
random (uniform distribution) permutation of {1,2,...,A^} (see |Ham72j ). Baik, Deift, 
and Johansson |BDJ99j show that 

(L(r) - 2r \ 
hm Pr M ; <s\ =F{s), 

where F(s) is the Tracy- Widom distribution for the largest eigenvalue of a GUE random 
matrix. The r^/^ in the denominator shows that x = 1/3 for this model, where x is the 
exponent describing the longitudinal fluctuations of a maximal increasing path. Moreover, 
Johansson jJohOObj applies the techniques from |NP95] and |LNP96] to show that the 
transversal fluctuations of the maximal paths have exponent ^ = 2/3, verifying the scaling 
identity x = 2^ — 1 for this model. 

3 Richardson's growth model and competition models 

3.1 The 1-type Richardson model 

First passage percolation can be put into the framework of interacting particle systems (see 
e.g. Lig85b| ) by defining a {0, 1}^ -valued process {r7i}i>o given by 



1 xiv^Bit), , , 

rit[v) = { ^ '' ioTve Z'^. 

I otherwise. 

We may think of sites in state as healthy and sites in state 1 as infected, so that the 
process represents an infection spreading outward from the origin. 
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When the passage times are i.i.d. exponentials with parameter A > 0, the memory- 
less property implies that the process r]t is Markovian. In this case the process is called 
Richardson's growth model [Ric73j , also known as the "contact process with no recover- 
ies," in comparison with the similarly defined contact process (cf. |Lig85b| ). In Richardson's 
growth model, a site in state 1 remains infected forever and tries to infect each of its 2d 
neighbors at rate A. Thus, the rate at which an uninfected site fiips from to 1 is equal to 
A times the number of infected neighbors it has. Only the origin is infected at time 0. 

As noted in |Ric73] . this process is related to a discrete time process called Eden's growth 
model |Ede61j . defined as follows: Set Ai = {0}, and forn > 1 set An = An-iU {vn} , where 
Vn is chosen from the set of uninfected sites with probability proportional to the number 
of neighbors it has in A„_i. Then An has the same distribution as B(tn), where 

tn = inf{t : B(t) contains n vertices}. 

3.2 The 2-type Richardson model 

In |HP98] . Haggstrom and Pemantle introduced the two-type Richardson model. In 
this model, instead of one type of particle spreading throughout the lattice, there are two 
species of particles competing for space. This competition is described by a {0, 1, 2}^ - 
valued Markov process {^t}t>o with parameters Ai and A2 which determine the fiip rates 
as follows: I's and 2's never fiip, while a fiips to a 1 (resp. a 2) at rate Ai (resp. A2) 
times the number of neighbors of type 1 (resp. 2). The I's and 2's represent sites infected 
by species 1 and 2, respectively, and O's represent uninfected sites. 

One natural question we can ask in the two-type model is whether both species continue 
growing indefinitely or whether one species ends up surrounded by the other so that it is 
only able to infect a finite number of sites. If Ai and A2 are two disjoint subsets of Z'^, we 
denote by Coex{Ai, A2) the event that both species eventually infect an infinite number of 
sites when species i initially occupies the sites in Ai, and we call this event coexistence 
or mutual unbounded growth for the initial configuration {Ai,A2). It is easy to see 
that Pr{Coex{Ai, A2)) < 1 unless both of the sets Ai and A2 are already infinite, so 
the first nontrivial question to ask is whether Pr(Coex(74i, ^42)) > 0. Clearly coexistence 
is impossible if one of the sets Ai surrounds the other set Aj, i.e. if there is no infinite 
path starting in Aj that does not intersect Ai. We say that the pair {Ai,A2) is fertile if 
neither set surrounds the other. Deijfen and Haggstrom showed that as long as the initial 
configuration of the process is finite and fertile, the choice of configuration is irrelevant to 
the question of whether coexistence has positive probability: 

Theorem 3.1 (Deijfen and Haggstrom |DH06aj ). If {Ai,A2) and {A[,A2) are two fertile 
pairs of disjoint finite sets in Z'^, then for any pair of growth rates Ai and A2, 

Pr{Coex{Ai, A2)) > ^ Pr{Coex{A[, A'^)) > 0. 

Weaker versions of Theorem 13.11 (in which the sets Ai and A^ consist of single points) 
appeared in |HP98j and [GM05J . and most treatments of coexistence have simply focused 
on the case where the initial configuration is ({0}, {I})- 

Intuitively, if the growth rates Ai and A2 are equal, we might expect that coexistence 
occurs with positive probability since neither species has an inherent advantage over the 
other. This result was in fact proved for d = 2 in |HP98j . and subsequently generalized to 
any d > 2 in |GM05j and |Hof05aj . On the other hand, if the growth rates are different, 
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say Ai > A2, then unless species 2 gets lucky and surrounds species 1 relatively quickly, 
species 1 is likely to overtake species 2 by virtue of its superior speed, making coexistence 
implausible. Haggstrom and Pemantle conjecture in |HP98] that Pr(Coex(0, 1)) = when 
Ai 7^ A2 and prove a somewhat weakened version of this conjecture in |HPOO] . In the next 
two sections we discuss the two-type Richardson model in more detail in the two cases 
Ai = A2 and Ai 7^ A2. 

3.3 Competition with equal growth rates 

The two-type Richardson model is somewhat simpler to analyze when both species grow at 
the same rate. If Ai = A2 = A, we can obtain the two-type process ^t as the projection of a 
single first passage percolation process with i.i.d. exponential(A) passage times, analogous 
to the definition of rjt in the one-type model. If the process starts with initial configuration 
(^1,^2), then 

{1 iiT{Ai,v) <t andT{Ai,v) <T{A2,v) 
2 iiT{A2,v)<t3.ndT{A2,v)<T{Ai,v) 
otherwise. 

Note that the definition of C,t can be generalized in an obvious way to model competition 
between k species with equal growth rates and initial configuration {Ai, . . . ,Ak), for any 
k > 1. Since the two-type (or fc-type) model and the one- type model are both defined in 
terms of an underlying first passage percolation process, results about one model can often 
be translated into results about the other, as will be illustrated below. 

Haggstrom and Pemantle first addressed the question of coexistence for species with 
equal growth rates in |HP98] . where they proved that Pr[Coex(0, 1)] > when d = 2. The 
main step in their proof was to show that in the related one-type process starting at 0, 
there are infinitely many sites in the right half-plane which have a > 50% probability of 
being infected after their neighbor to the left, so that these sites "sense" that the infection 
is coming from the left. From there, it is a small step to show that in the two- type process, 
with positive probability there are infinietly many sites in the right half plane that are 
reached by species 2 before they are reached by species 1, and that a symmetric situation 
holds in the left half-plane. 

Observe that the definition of C,t makes sense for more general passage time distributions, 
although the Markov property holds only in the i.i.d. exponential case. However, for any 
stationary distribution of passage times, the two species will still be growing at the same 
average rate, and we might expect coexistence to hold in the stationary case as well. 
Indeed, Caret and Marchand |GM05j and Hoffman |Hof05a] independently generalized the 
coexistence result of |HP98j to a large class of stationary ergodic passage times in any 
dimension d > 2. Furthermore, an analogue of Theorem 13. II holds in the stationary case so 
that the starting configuration is still irrelevant |GM05t p. 312]. 

Coexistence in the two-type or fc-type model is related to the existence of one-sided 
geodesies in the corresponding one-type model. If coexistence of k species occurs, then 
the same compactness argument used to show that K{r{0)) > 1 shows that there exist k 
disjoint one-sided geodesies in the underlying first passage percolation process, one starting 
in each of the initial sets Ai, . . . ,Ak. Therefore, denoting coexistence in the /c-type model 
by Coex(y4i, . . . , Ak) and the existence of disjoint geodesies Gi starting in the sets Ai by 
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Geo(y4i, . . . , Ak), we have 

Pr[Geo(A, . . . , Ak)] > Pr[Coex(A, . . . , A^)]. (3.3.1) 

Furthermore, if Geo{Ai, . . . , Ak) occurs, it seems plausible that some finite modification 
of passage times might allow the construction of k one-sided geodesies starting at so 
that Pr[i^(r(0)) > A;] > 0. In fact, at least when k = 2 (and probably for any k - see 
[HofOSbj ) . we can go in the other direction as well, from geodesies to coexistence: For the 
class of stationary measures considered in |GM05] or |Hof05aj . it can be shown (see [GM05 
Lemma 5.3]) that 

Pr[Coex(0, 1)] > ^ Pr[K{T{6)) > 2] > 0. 

Thus, since coexistence of two species occurs with positive probability, there are at least 
two one-sided geodesies starting at with positive probability. In fact, while Garet and 
Marchand |GM05] use techniques similar to those in |HP98j to prove that coexistence is 



possible and then conclude that there are at least two geodesies with positive probability, 
Hoffman |Hof05aj first proves that there almost surely exist at least two distinct one- 
sided geodesies (not necessarily with the same starting point) and uses this to show that 
coexistence has positive probability. 

In |Hof05bj . Hoffman applies the techniques in |Hof05a] to the fc-type model to obtain 
further results about both geodesies and coexistence when d = 2 and the passage times are 
given by a certain class of "good" ergodic stationary measures z/. Although the results are 
stated only for d = 2, the methods can be applied to any d > 2. We now state the main 
results, which depend on the geometry of the limit shape Bq corresponding to u. For a 
good measure z/, let Sides(i/) be the number of sides of dBo if dBo is a polygon or infinity 
if dBo is not a polygon. 

Theorem 3.2 (Hoffman |Hof05bj ). Let u be a good stationary measure on (]R+)^*-^ ', and 

let k < Sides(z/). For any e > 0, if r is sufficiently large there exist ui, . . . ,Uk € d{rBo) 
such that 

Pr[Coex(Mi, . . . , Uk)] > 1 — e and (by (13.3. ip ) Pr[Geo(ui, . . . , Uk)] > 1 — e. 

The points Ui in Theorem 13.21 are chosen to be the lattice points closest to points 
u[, . . . ,u'i. G dirBo) at which the tangent fines of dirBo) are distinct (such points exist by 
the assumption that k < Sides(i/)). The idea of the proof is that if f i, . . . , Vk are points 
on OBq with distinct tangent lines L^^, . . . , L^^., then with positive probability, for each i 
there will be infinitely many n such that Vi is closer (in travel time) to the translated line 
nvi + Ly- than any of the other points Vj are. This shows that if the process starts with 
initial configuration {vi, . . . , Vk) (assuming Vi G Z^), each Vi will infect infinitely many sites 
with positive probability, so coexistence occurs. By scaling the picture up by a sufficiently 
large factor r, the probability of coexistence can be made arbitrarily close to 1. 

Hoffman also obtains the following results about one-sided geodesies starting at 0. 

Theorem 3.3 (Hoffman |Hof05bj ). Let u be a good stationary measure on (R^)^^'^ \ If 

k < Sides(z/)/2, then 

K{T{6)) > k a.s. 

Theorem 3.4 (Hoffman |Hof05bj ) . Let v = (C{t))'^^^'^ \ where r is an exponential random 
variable. If k < Sides(i/), then 

Fr[K{T{0)) >k]>0. 
17 



Observe that by symmetry, we must have Sides (z/) > 4 when li = 2, so Theorem 13.21 
shows that coexistence of four species is possible for any good measure u, and Theorem 13.41 
shows that with i.i.d. exponential passage times it is possible to get four one-sided geodesies 
starting at 0. Furthermore, in |DL81] it is shown that there is a nontrivial i.i.d. measure u 
such that Bq is neither a square nor a diamond, so by symmetry we must have Sides(z/) > 8. 
Thus, Theorem 13.21 implies that there is a nontrivial i.i.d. measure z/ for which coexistence 
of eight species is possible. In |HM95j it is shown that there exists a good measure u such 
that Bq is the unit disc, so Theorem 13.31 implies that there exists a good measure z/ such 
that K(T{6)) = oo a.s. 

3.4 Competition with different growth rates 

If the growth rates Ai and A2 are different for the two species, we can construct the Markov 
process ^t from two independent i.i.d. exponential first passage percolation processes on Z'^, 
one with parameter Ai and the other with parameter A2. However, the description of C,t is 
not quite as simple as it was in the case of equal growth rates because, with two underlying 
sets of passage times instead of just one, there is no guarantee that the geodesies between 
infected vertices of one species will not cross geodesies of the other species. For this reason, 
the values of ^t must be defined iteratively by considering the process at the time tn of 
the n*'^ infection. We mention that the state of the process at time tn can be described 
analogously to Eden's growth model for the one-type process, except that the edge causing 
the next infection is chosen from all edges on the boundary with probability proportional 
to Aj if the edge borders species i. 

In |HP98] . Haggstrom and Pemantle conjecture that coexistence in the two-type Richard- 
son model is impossible if Ai 7^ A2. While the full conjecture is still an open problem, 
Haggstrom and Pemantle were able to prove the slightly weaker result that coexistence is 
impossible for almost all choices of parameter values: 

Theorem 3.5 (Haggstrom and Pemantle |HPOO] ). For the two-type Richardson model on 
Z"^, d > 2, with Xi = 1 we have 

Pr(Coex(0, 1)) = 

for all but at most countably many choices of \2. 

By time scaling, the probability of coexistence depends only on the ratio A = A2/A1, so 
Theorem 13.51 remains true for any other choice of Ai. Furthermore, the same time-scaling 
argument plus symmetry implies that the probabilities of coexistence for the pairs (1,A) 
and (1, 1/A) are equal, so it suffices to consider the case Ai = 1 and A2 = A G [0, 1]. 

At first glance, it may seem strange that Theorem 13.51 has not been extended to in- 
clude all values of A2 7^ Ai. Intuitively, we expect that Pr(Coex(0, 1)) should decrease as 
A = A2/A1 moves farther away from 1. Since Theorem 13.51 implies that we can choose A 
arbitrarily close to 1 such that Pr(Coex(0, 1)) = 0, such monotonicity would imply that 
coexistence is impossible for all A 7^ 1. However, it is not obvious how to prove that the 
probability of coexistence is monotone in A. In fact, although this monotonicity property is 
certainly plausible for the integer lattice 7/, Deijfen and Haggstrom |DH06b] have shown 
that there are other (highly non-symmetric) graphs where monotonicity does not hold. 

We now give a brief outline of the proof of Theorem 13.51 from |HP00j . The main tool 
is the following proposition, which we state as it appears in |DH06cj . Let Y* ^^ ^^ denote 



the law of the two-type process with rates Ai, A2 and initial configuration (^11,^2). For 
i = 1,2, let Gi be the event that species i finally infects an infinite number of sites (so 
Coex(y4i, y42) = Gi n G2), and let Bq denote the limit shape for the one-type Richardson 
model with rate 1. Then 

Proposition 3.6 f |HP00l Prop. 2.2], |DH06cl Prop. 5.2]). For any X < 1 and e > we 

have 

lim sup P^^Aal'^s) = 0, 

where the supremum is over all initial configurations {Ai,A2) such that 

A2 is contained in tBq, while 

Ai is not contained in (1 -|- ejrB^. (3.4.1) 

For example. Proposition 13. 61 says that if we start the process with with the slow species 
occupying the entire /i-ball of radius r and the fast species occupying a single site outside 
the /i-ball of radius (1 + e)r (where /i is the norm for the unit rate Richardson model), the 
survival probability of the slow species goes to zero as r — )■ 00. Using Proposition 13.61 and 
the strong Markov property, Haggstrom and Pemantle show that if Coex(0, 1) occurs, the 
set of sites infected by both species, scaled by t, converges a.s. to the limit shape of the 
slow species. 

To prove Theorem 13. 5^ Haggstrom and Pemantle first describe a coupling Q of the 
processes P^'j. (A G [0, 1]) such that Q-a.s., for all t, the set I's at time t decreases with A, 

and the set of 2's at time t increases with A. Writing Coex(A) for the event that Coex(0, 1) 
occurs at parameter A under the law Q, they use the above result to show that Q-a.s., 
Coex(A) occurs for at most one A G [0, 1]. That is, 

Q (lcocx(A) = for all but at most one A G [0, 1]) = 1, 

so by Fubini's theorem, 

Y, Q(Coex(A)) = Eq 5^ lcoex(A) < Eq 1 = 1. 

ag[o,i] Ae[o,i] 

Therefore, since the sum on the left is finite, there can be only countably many A's such 
that Q(Coex(A)) > 0, which proves Theorem 13.51 since Q(Coex(A)) = P^':j.(Coex(0, 1)). 

While the results in |HPOO] apply only to finite initial configurations, Deijfen and 
Haggstrom [DH06cj recently proved some interesting results about coexistence in the case 
when one of the initial sets Ai is infinite. In particular, they considered the cases where 
Ai is either the hyperplane H = Hq = {z & TJ^ : Zi = Q} (minus the origin) or the half 
line L = {z ^U^ : Zi <Q and Zi = for alH 7^ 1} (minus the origin), and A2 = {0}. Their 
main result is 

Theorem 3.7 (Deijfen and Haggstrom |DH06cj ). For the two-type Richardson model in 
d >2 dimensions, 

1. Pr[Coex(/7 \ {0}, 0)] > z/ and only if Ai < Aa- 

2. Pr[Coex(L \ {0}, 0)] > «/ and only if Ai < Aa- 
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The fact that coexistence is impossible if Ai > A2 for either H 01 L follows from 
Proposition 13.61 since whenever Ai is infinite while A2 is finite, the pair {Ai,A2) satisfies 
(13.4. ip for all sufficiently large r. The "if direction of[T]is proved by combining a shape 
theorem for the one-type process starting from H (which is proved using the large deviation 
bounds of Theorem l2.6p with a shape theorem for the "hampered" one-type process starting 
from and restricted to a cylinder about the first coordinate axis (which follows from a 
standard modification of the proof of the ordinary shape theorem). The strategy of proof is 
to show that when Ai < A2, there is a positive probability that species 2 gets a big enough 
head start over species 1 that it is able to take over the entire cylinder without interference. 
The corresponding result for L follows from the result for H because (a rotation of) L is a 
subset of H, and the probability of survival for either species is monotone with respect to 
the starting configuration [DHOGcl Lemma 3.1]. 

For the critical case Ai = A2, the proof that coexistence is possible when Ai = L\ {0} 
follows techniques similar to those used in |HP98] . In fact, the coexistence result for 
Ai = L\ {0} when Ai = A2 allows an easy proof of the the coexistence result in [HP98] 
(see [DHOGcj Theorem 6.1]). The proof that coexistence is impossible when Ai = H \ {0} 
is rather more involved, and we will not discuss it here. We mention, however, that this 
result shows that in a one-type process started from H, almost surely every vertex in H 
will infect only a finite number of vertices in Z*^. 

As was the case with equal growth rates, the definition of the process ^t in terms of first 
passage percolation makes sense for more general passage times, although again, Markovity 
will be lost in the non-exponential case. In |GM06j , Garet and Marchand extend the results 
of |HP00j to include i.i.d. passage times which are not necessarily exponential but for which 
the passage time distributions for the two species are stochastically comparable. In this 
setting, they show that for any d, if the slow species survives, the fast species cannot occupy 
a very high density of space (for example, "it could not be observed by a medium resolution 
satellite"). For d = 2, they show that almost surely, one species must finally occupy a set 
of full density in the plane while the other species occupies only a set of null density. They 
also obtain deviation bounds similar to those in [HPOO] showing that if coexistence occurs 
then the infected region in the two-type process must grow essentially according to the law 
of the first passage percolation process governing the slow species. Finally, they prove an 
analogue of Theorem 13.51 for families of stochastically comparable passage times indexed 
by a continuous parameter, showing that coexistence cannot occur except perhaps for a 
countable set of parameters. 
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